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Abstract. The statistics of the nodal lines and nodal domains of the eigenfunctions 
\q | of quantum billiards have recently been observed to be fingerprints of the chaoticity of 

the underlying classical motion by Blum et al (2002 Phys. Rev. Lett. 88 114101) and 
by Bogomolny and Schmit (2002 Phys. Rev. Lett. 88 114102). These statistics were 
shown to be computable from the random wave model of the eigenfunctions. We here 
study the analogous problem for chaotic maps whose phase space is the two-torus. We 
show that the distributions of the numbers of nodal points and nodal domains of the 
eigenvectors of the corresponding quantum maps can be computed straightforwardly 
and exactly using random matrix theory. We compare the predictions with the results 
of numerical computations involving quantum perturbed cat maps. 
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In a recent article Blum et al (2002) observed that the number-distributions of the nodal 
domains of quantum wavefunctions of billiards whose classical dynamics is integrable 
are different from those for chaotic billiards and argued that the latter are universal. 
Thus, the number-distribution of nodal domains appears to be a new criterion for 
^ ■ quantum chaos that complements the usual ones based on spectral fluctuations. Blum 
et al computed these distributions for some integrable (and separable) systems, but no 
analytic formula exists for the number of nodal domains of a chaotic billiard. Berry 
(1977) has conjectured that the wavefunctions of quantum systems with a chaotic 
classical limit behave like Gaussian random functions. Supported by numerical evidence, 
Blum et al found that the limiting distribution of the number of nodal domains can be 
reproduced assuming Berry's conjecture. Bogomolny and Schmit (2002) developed a 
percolation model for nodal domains of Gaussian random functions and showed that 
their number is Gaussian distributed. They computed the mean and variance of this 
distribution, which are both proportional to the mean spectral counting function. Their 
results agree with the numerical computations reported by Blum et al for chaotic 
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billiards. The influence of a boundary on the nodal lines of Gaussian random functions 
has been investigated by Berry (2002), Gnutzmann et al (2002), and Berry and Ishio 
(2002). This is expected to model the nodal properties of billiard wavefunctions near 
boundaries. 

We here consider the analogous problem for one- dimensional time-reversal- 
symmetric systems with discrete time evolution and whose phase space is the two- 
dimensional torus T 2 . The classical dynamics of such systems corresponds to the 
action of symplectic maps on T 2 , and their quantum mechanics to that of unitary 
matrices Un (called propagators or quantum maps) on a Hilbert space of dimension 
N — 1/h, where h is Planck's constant. Modelling the eigenvectors of Un by those of 
random unitary symmetric matrices (such matrices constitute the circular orthogonal 
ensemble, COE, of random matrix theory), we compute the number-distributions of 
nodal domains and nodal points (the analogues of nodal lines in billiards) exactly. It 
is shown that these become Gaussian as iV — > oo and that the mean and variance are 
proportional to iV (precisely as in the billiard case). We compare our results with 
numerical computations involving the eigenvectors of perturbations of quantum cat 
maps whose classical dynamics are hyperbolic and whose spectral statistics are known 
to be accurately predicted by random matrix theory (Basilio de Matos and Ozorio de 
Almeida 1995, Keating and Mezzadri 2000). 

Consider the Helmholtz equation with Dirichlet boundary conditions 

- A#(r) = EV(t), rGfl, (1) 

where Q is a connected compact domain in a two-dimensional Riemann manifold. The 
nodal lines are the zero sets of real solutions of equation (JTJ); the nodal domains are 
connected domains in Q where \l/(r) has constant sign. Now, let {^n( r )K£Li De a set 
of eigenfunctions of the laplacian on Q ordered by the magnitude of the corresponding 
eigenvalue E n , and let u n be the number of nodal domains of the n-th eigenfunction. 
Courant (1923) proved that v n < n. Let I g (E) = [E, E + gE], for g > 0. Blum et al 
(2002) introduced the distribution 

P,{xJ 9 {E)) = j r £ 5 { x -^)> ( 2 ) 

1 E n el g (E) 

where iVj is the number of energy levels in I g (E). The limiting distribution of nodal 
domains is defined by 

P h (x) = lim P h (x,I g (E)). (3) 

We now introduce a density that is the analogue of (JHJ) for quantum maps. The 
periodicity of the two-torus constrains the wavefunction to be an infinite sum of delta- 
functions supported at rational points of the form j/N, with j integer, in both the 
position and momentum basis (Hannay and Berry 1980), i.e. 

N ( ' \ 
^(l) = EE C ^ ( q ~ a7 + m ) ' ( 4a ) 
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N , 



m£Z j=l 

where N — 1/h and 



1 / i 

Moreover, since VK?) an d ^(p) are periodic, Cj = Cj + jy and Cj = Cj + jy. Therefore, a 
quantum state is completely determined by N complex numbers, which implies that the 
Hilbert space is isomorphic to C^. The coefficient Cj can thus be interpreted as the 
value of ip(q) at q — j/N (the Heisenberg uncertainty principle is not violated, because 
the periodic sum of delta-functions that defines ip{p) extends to infinity). Now, let Un 
be the matrix realization of a quantum map in the basis {\j )}f=i, where 



g\j) = 



J2 5 (^-i+ m )- ( 6 ) 



We shall consider only systems whose dynamics is invariant under time reversal, so that 
Un is a symmetric unitary matrix and, without loss of generality, the eigenvectors can 
be taken to be real. 

Because of the topology of the phase space, an eigenvector of C/jv is equivalent to a 
sequence of N real numbers with periodic boundary conditions, i.e. c\ = cn+i- A nodal 
point is then identified whenever two consecutive coefficients Cj have opposite sign. The 
total number of nodal points in a given eigenvector is 

1 N 

v = 2 - s S n ( c i) s g n ( c i+i)]- ( 7 ) 
i=i 

where 

{1 if x > 0, 

if x = 0, (8) 

-1 ifx<0. 

Similarly, a nodal domain is a set of consecutive integers + j + 2, . . . , j + k} such that 
the corresponding coefficients Cj lie between two nodal points and thus have constant 
sign. As a consequence of the periodicity of the coefficients Cj, there can be only an even 
number of nodal points, equal to the number of nodal domains; the only exception is 
when there are no nodal points and only one nodal domain. It follows from the results 
to be presented later that as N — > oo the probability that all the Cj s have the same 
sign is negligible, and so we shall denote by v both the number of nodal points and the 
number of nodal domains. Finally, the limiting distribution is defined by 

N 



n=l 

where, as for billiards, v n is the number of nodal domains (points) of the n-th eigenvector. 
Identical definitions can obviously be formulated in the momentum representation. 
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When the classical limit of Un is a chaotic map, the eigenstatistics of Un are 
expected to be the same as those of matrices in the COE (Bohigas et al 1984). The 
COE probability measure is invariant under the mapping 

U i-> OUO T , (10) 

where U is a unitary symmetric matrix and O is an arbitrary orthogonal matrix. Hence, 
each eigenvector of U is mapped by an orthogonal transformation into an eigenvector 
of a new matrix that by (fTU|) has the same weight in the ensemble as U and the same 
spectrum. As a consequence (see, e.g. Haake 2000), the eigenvectors of matrices in 
the COE are uniformly distributed on the unit sphere in Wt N and the joint probability 
density of their components is 

P OOE ( Cl , C2 ,..., CAr ) = -A-r^V( i-> "A ) • «id 




v / v j= i 

The above distribution is independent of the signs of the Cj s, therefore they can be either 
positive or negative with equal probability and there are no correlations among the signs 
of different coefficients. This simple observation allows us to compute analytically all 
the relevant quantities in a very straightforward way. 

The signs of the Cj s behave like a sequence of N independent random variables 
Sj that assume the values {1,-1} with equal probability ~; in other words, they are 
equivalent to an array of non-interacting particles with spin | and periodic boundary 
conditions. Thus, the probability of a configuration with N + spins up and AL = N — N+ 
spins down is given by the binomial distribution 

P(N + ,N-.) = ±( J V (12) 




The computation of the density (J9J) requires a simple combinatorial argument. In a 
periodic chain of N spins there are N possible positions where a nodal point can be 
located. Hence, the number of configurations with v nodal points is zero when v is 
odd and twice the number of ways of choosing v objects among N, irrespective of their 
ordering, for even u, i.e. 

[1 + (-1)1 V (13) 

The factor of two in front of the binomial coefficient is due to the fact that by 
simultaneously changing the sign of all the spins in the chain we obtain a new 
configuration with the nodal points in the same positions. Finally, the distribution 
of the number of nodal points and nodal domains is given by 

W) =WL{» ]. (14 ) 



V 



The mean (y) and variance a = (y ) — (v) can be easily computed: 

("> = y and a 2 = ^. (15) 
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Figure 1. Nodal domain distributions (0) of the unperturbed quantum cat map (a) 

and of its perturbation (|20|l with k = 0.30 (b) compared with the Gaussian l|16f) ( )■ 

The dimension of the Hilbert space is N — 3511. 



Equations ()14j) and (fl3|) correspond to the results that Bogomolny and Schmit (2002) 
obtained for the percolation model of random wave functions in two-dimensional 
systems. By letting iV — ► oo and scaling x = u/N, the discrete distribution (JT4*j) tends 
to a continuous Gaussian probability density with mean | and variance a 2 = l/4iV, i.e. 



P m (x, N) ~ y~~ ex P [- 2N ( x ~ V 2 ) 2 ] , iV ^ oo. (16) 

This is the main result of this note. 

In order to compare the distribution (fT7)|) with numerical computations, we consider 
perturbations of the following hyperbolic (cat) map: 



A : \ q ~ I I M mod 1. (17) 




Because of the number-theoretical properties of A, the spectrum of the propagator 
Un(A) is non-generic (Keating 1991, Kurlberg and Rudnick 2000) in that it does not 
obey the random matrix theory conjecture. However, if a small nonlinear perturbation 
is introduced, the composite map is still hyperbolic but loses its arithmetical nature. As 
a consequence, the spectrum of the new quantum map has random matrix correlations. 
Hence, we perturb fll7j) with the following shear in the momentum 

P: (p)"(p+£»s(2^)) <18) 
and study the propagator Un(4>) of the map 

<p = poAop. (19) 
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U N (4>)i m = -= exp 
ViiV I 



^ - Im + + -^5- (sin(27r//iV) + sin(27rm/iV)) 



(20) 



The matrix elements of this propagator in the basis (JEJ) are 

2vri 

r — Im + + - 

(Basilio de Matos and Ozorio de Almeida 1995). It can be shown that the only symmetry 
of this quantum map is time reversal (Keating and Mezzadri 2000). Furthermore, 
if k < k max = 0.32..., then the map (fT9*|) is uniformly hyperbolic and the spectral 
statistics of the propagator (j20j) are consistent with random matrix theory (Basilio de 
Matos and Ozorio de Almeida 1995). Figure |T(b)| shows the nodal domain distribution 
of the eigenvectors of the quantum map ()20|) for a particular choice of k and N, together 
with the density (jlfij) . The nodal domain distribution of the unperturbed quantum map, 
figure [T(a)| also appears to be Gaussian, but its variance cannot be predicted by random 
matrix theory. 

As the perturbation parameter k varies, the nodal points in a given eigenvector 
of the matrix (|20p change their positions. A natural question then arises: what is the 
minimum number of parameters needed to create or coalesce nodal points and alter the 
number of nodal domains? In other words, what is the codimension of the nodal points? 
Since the spins in a chain are uncorrelated, the functions Sj(k) will be independent, 
therefore nodal points move randomly without repelling or attracting each other. Thus, 
the codimension of nodal points is one and a single parameter is enough to create or 
annihilate nodal domains with equal probability. This behaviour is illustrated in figureEl 
the scaled number of nodal domains x(k) = v{k)/N of an eigenvector oscillates around 
|, and since the Sj(k) are independent, the value distribution of x{k) is given by the 
Gaussian (|TB|) . 

Finally, it is worth remarking that this problem is equivalent to a one- dimensional 
Ising model of non-interacting spins in a magnetic field B with periodic boundary 
conditions, whose Hamiltonian and partition function are 

N 

H = -B^2 s i> s j = ±l ' s i = s n+i (21) 

and 

Z(P, B) = ]T ■ ■ E ex P (-W = 2N ™MPB) N (22) 

{si} {s 2 } {s N } 

respectively. All the relevant thermodynamical quantities should be computed at 
/3 — B — 0. This plays the role in this case of the analogy between the nodal statistics of 
billiard wavefunctions and the Potts model suggested by Bogomolny and Schmit (2002). 
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Figure 2. (a) Scaled number of nodal domains x(k) — v(k)/N of an eigenvector 
of the matrix J3UJ), with N = 1069, as a function of the perturbation parameter k; 
(b) value distribution of x(k) averaged over all eigenvectors (0) compared with the 
Gaussian (fTB)l ( ). 
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